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We show that quantum correlations induced by electron-electron interactions in the presence of 
random impurity scattering can play an important role in the thermal stabilization of amorphous 
Hume-Rothery systems: When there is strong backscattering off local, concentrical ion clusters, the 
static electron density response x(0, q) acquires a powerlaw divergence at q = 2&f even at elevated 
temperature. This leads to an enhancement as well as to a systematical phase shift of the Friedel 
oscillations, both consistent with experiments. The possible importance of this effect in icosahedral 
quasicrystals is discussed. 

I. INTRODUCTION II. THEORY 



A large class of noble-polyvalent metal alloys exhibit 
a crystalline to amorphous transformation (CAT) as a 
function of the polyvalent metal content. There are sev- 
eral experimental indications for the amorphous phase 
being stabilized by the Hume- Rothery (HR) mechanism, 
i.e. by forming a structure-induced pseudogap in the 
electronic density of states (DOS) at the Fermi level ep. 
The experimental evidence jjj includes the observation 
of a pronounced pseudogap, a maximum of the electrical 
resistivity at the CAT, and the coincidence of the atomic 
spacing a with the Friedel wave length Xp = ir/kp near 
the CAT, where hkp is the Fermi momentum. Very simi- 
lar behavior is found in icosahedral (i) quasicrystals 
The conjecture of a HR-like stabilization mechanism is 
strongly supported by detailed theoretical studies both 
for amorphous Q and for quasicrystalline J|-|§ systems, 
although there is also the possibility of entropic stabiliza- 
tion H . The fact that in the amorphous state the above- 
mentioned structural matching is observed over distances 
of up to 5Af has lead to the assumption that these re- 
markably long-range correlations are induced by the ions 
being bound in the minima of the potential formed by the 
Friedel oscillations (FO) around an arbitrary central ion. 

However, several experimental puzzles have remained 
unexplained: (1) At finite temperature T and also in the 
presence of disorder the impurity-averaged FOs are expo- 
nentially damped due to the spread of the Fermi momen- 
tum over a width given by T and the inverse elastic mean 
free path, respectively. Hence, the stability of amorphous 
alloys at elevated T, in particular their intermediate- 
range structural correlations, are difficult to explain by 
the conventional FOs. (2) In all HR-systems, invariably 
the amorphous state is thermally most stable just at the 
CAT Q|. (3) In the amorphous state the ionic positions 
are systematically shifted compared to the minima of the 
Friedel potential of a free electron sea [0. The experi- 
mental findings (l)-(3) raise the question of a system- 
atical, composition dependent enhancement and phase 
shift if of the FOs, where ip varies from tp ~ 7r deep in- 
side the amorphous phase to ip — ir/2 at the CAT. As 
shown below, these problems can be explained by one 
single quantum effect. 



In amorphous alloys the electronic motion is diffusive 
instead of ballistic. Since diffusion, as a dissipative pro- 
cess, is difficult to include in an ab initio calculation, we 
here choose a Feynman diagram technique, where diffu- 
sion arises in the formalism by averaging over all random 
configurations of the system. The Nagel-Tauc condition 
2k p ~ 2n/a = k p implies strong electronic backscatter- 
ing off local ion clusters. It has been shown JTo| [l2|| that 
this not only generates a pseudogap but at the same time 
leads to a substantial enhancement of the electron trans- 
port or density relaxation rate r _1 over the quasiparticle 
decay rate . This is evidenced experimentally by the 
anomalously small electrical conductivity, a = ne 2 r/m*. 
Thus, we have as a generic feature of the amorphous 
state: r _1 » t~^. 

We now turn to the calculation of the electron density 
distribution p(r) around an ion embedded in the disor- 
dered electron sea. One does not expect a single ion to 
generate a sizable phase shift of the FOs at distances 
r > a, p(r) oc cos(2kpr — f)/r 3 , since it is equal to the 
electron backscattering phase off that ion, which is small 
unless £p is close to an internal resonance. Instead we 
consider quantum effects due to disorder and Coulomb 
interaction and assume, for simplicity, a point-like ion 
charge and a quadratic band Sk — h 2 k 2 /2m*. In a diffu- 
sive electron system screening is inhibited, so that the ef- 
fective Coulomb interaction <o^* (z, Z) between electrons 
with complex frequencies z and z + Z acquires a long- 
range, retarded part Q , 

vff(z,Z) = - wp ^-^ T 2 (z,Z,q), «„ = ^p-, (1) 

where e RPA (Z, q) = l+2Tticj j{Z sgaZ" +iq 2 D) is the dis- 
ordered RPA dynamical dielectric function and the dif- 
fusion vertex, defined in Fig. [l] a), is 



l/T S ^ Z " z'\z + zy < o 



T(z, Z,q) = < z+iq*D sgnz 

1 otherwise 



(2) 



D = 1/3 vpr and " denote the diffusion constant and 
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FIG. 1. a) Diffusion vertex I\ b) Polarisation 11(0, q) in- 
cluding leading order quantum correction induced by disorder 
and interactions. Dashed lines denote electron-ion scattering, 
the wavy line with solid triangles the effective Coulomb inter- 
action. 

the imaginary part, respectively. The long-range nature 
of Vqf* originates from the diffusion pole of T. In or- 
der to calculate its effect on the FOs, one must consider 
contributions to the polarisation 11(0, q) where T(z, Z, q') 
enters at Z, q 1 ~ 0, although 11(0, q) is evaluated at large 
external momenta q ~ 2k F . The leading singular con- 
tribution arises from the quantum correction shown in 
the second diagram of Fig. [j] b) . For amorphous metals 
(t -1 ^> Tqp 1 , see above) it is evaluated as [ fL2| , 



nW(o,?) = C( E n) / 



i) 



cosh ^ ^J\x — 1| 



(3) 



where x — x(v) = (q/2k F )/ y/l + vje F and C(e f t) — 
-0.343 [2m*k F /(2nh) 2 }/(e F T) 7 / 2 . The first term of 
Fig. |l| b), n(°)(0, g), corresponds to the Lindhard func- 
tion, where T contributes only a nonsingular factor of 
0(1). It is seen from Eq. (§) that for T = 0, n< 1 )(0,g) 
exhibits a powerlaw divergence oc — sgn(g — 2k F )/\q — 
2k F \ 1 ' 2 at q = 2k F . Although at finite T or r^ 1 the 

divergence of EIW(0, q) is reduced to a peak, the in- 
verse dielectric function l/e(q) = 1/(1 — v q ll(0,q)) still 
has a q — 2k F divergence at a critical transport rate 
t~ (T) even for non-zero single-particle relaxation rate 
T qp < T_1 an d &t finite T. The parameter r _1 is varied 
experimentally by changing the composition of the alloy. 



III. DISCUSSION AND COMPARISON WITH 
EXPERIMENTS 



Fourier transforming 1 — l/e{q) to obtain p(r) fllCf ] 
shows that for incomplete Fermi surface-Jones zone 
matching, i.e. small r _1 , the quantum corrections gener- 
ate density oscillations p^(r) oc — cos(2fci?r)/r 3 , which 
overcompensate the conventional FOs, implying a phase 
shift of ip = 7r [Q. As r _1 — * t^ 1 -, the increasing 2k F 
peak of l/e(q) leads in addition to density oscillations 
p^\r) oc sin(2k F r)/r 2 , so that in the vicinity of r^ 1 



p(r) 



cos(2k F r) 
{2k F rf 



sm(2k F r) 
' {2k F r) 2 ' 



with A^- 1 ) ~ 0.343tt(1 - r" 1 /^ 1 )- 1 / 2 . The expo- 
nent 1/2 is characteristic for diffusive behavior. Thus, 
the FOs are shifted by ip = ir — tan [2fe.Fr A(r -1 )] ~ 
tt/2 + l/(2k F rA), i.e. the diverging Friedel amplitude 
necessarily goes hand in hand with ip = tt/2. Note that, 
in contrast to the conventional FOs, this divergence is 
robust against damping due to finite T or disorder. The 
point where the amplitude A diverges should be identi- 
fied with the CAT, since at this point the fluctuations 
of the Friedel potential also become large, allowing the 
system to find its crystalline ground state. This resolves 
in a natural way the problems (l)-(3) mentioned in the 
introduction. 
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FIG. 2. Crystallization temperature T c as a function of 
the DOS suppression at ef, (1 — <?)■ Data points represent T c 
for a,-Cu\- x Sn x [1]. The solid curve is the fit of the present 
theory (see text). Vertical line: position of CAT. The inset 
shows the phase shift p of the first maximum of the charge 
density distribution p(r). Solid line: theory. Data points with 
solid line: measurements [1] for &-Cu\- x Sn x . 



For a direct comparison with experiments the con- 
trol parameter of the theory, r _1 , must be translated 
into a parameter, which is experimentally accessible: It 
follows from the scattering theory |pd| , |l2| that t _1 = 
t^ 1 + 7 S(2k F ), where the peak of the ionic structure 
factor S(q = 2k F ) controls the backscattering amplitude, 
7 is a constant, and t^ 1 is an offset due to momentum 
independent scattering. S(2k F ) in turn is proportional 
jjj to the measured, structure-induced suppression of the 
DOS N(e F ) at the Fermi level, l~N(e F )/N (e F ) = 1-g, 
compared to the free electron gas, N (s F ). The resulting 
fit of the crystallization temperature T c is shown in Fig. 

where the contribution to the stability coming from the 
pseudogap formation is assumed to be linear in (1 — g) 
(dashed line). The inset shows the calculated phase shift 
(p and the measured shift of the atomic nearest neighbor 
position relative to the position of the first conventional 
Friedel minimum, a — ir/kp. Note that there is no ad- 
justable parameter in ip. The general behavior of the shift 
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is well explained by the theory; however, the experimen- 
tal data approaches cp = ir/2 faster than predicted. This 
might be attributed to the fact that, as seen from the 
discussion after Eq. (Q), the higher-order Friedel min- 
ima approach <p = it/2 faster than the first one. In this 
light, the agreement between theory and experiment is 
remarkably good. 

The structural similarities H between amorphous al- 
loys and «-quasicrystals suggest that the quantum effect 
discussed above may be important in the latter systems 
as well. In fact, quasicrystals seem to fulfill all the nec- 
essary preconditions for this effect to occur, i.e. effec- 
tively diffusive electron motion [|l4| and r _1 » t^ 1 . 
The latter is supported by the Fermi surface matching, 
i.e. by the experimental observation || and theoretical 
prediction of structure-induced pseudogaps. More- 
over, another more commonly known effect of disorder- 
enhanced Coulomb interaction, the y/\E — Ef\ behavior 
of the DOS in the pseudogap [|l3|, may have been already 
observed in i-quasicrystals by tunneling measurements 
of the DOS ^j. It is proposed to include the enhanced 
Friedel potential calculated in the present work in the 
pseudopotential of more quantitative ab initio calcula- 
tions. 
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